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ABSTRACT 

NGC 7217 is an unbarred early- type spiral galaxy having a multi-segment exponential 
light profile and a system of starforming rings of the unknown origin; it also possesses 
a circumnuclear gaseous polar disc. We analysed new long slit spectroscopic data for 
NGC 7217 and derived the radial distributions of its stellar population parameters 
and stellar and gaseous kinematics up to the radius of r « 100" (~ 8 kpc). We 
performed the dynamical analysis of the galaxy by recovering its velocity ellipsoid at 
different radii, and estimated the scalehcights of its two exponential discs. The inner 
exponential stellar disc of NGC 7217 appears to be thin and harbours intermediate 
age stars (<ssp ~ 5 Gyr). The outer stellar disc seen between the radii of 4 and 7 kpc 
is very thick (z = 1 . . . 3 kpc), metal-poor, [Fe/H]< —0.4 dex, and has predominantly 
young stars, tssp = 2 Gyr. The remnants of minor mergers of gas- rich satellites 
with an early-type giant disc galaxy available in the GalMer database well resemble 
different structural components of NGC 7217, suggesting two minor merger events in 
the past responsible for the formation of the inner polar gaseous disc and large outer 
starforming ring. Another possibility to form the outer ring is the re-accretion of the 
tidal streams created by the first minor merger. 
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1 INTRODUCTION 

Currently, our views on galaxy structure are changing dra- 
matically. The classical view that every disc galaxy rep- 
resents a combination of a de-Vaucouleurs' bulge, a sort 
of an elliptical galaxy inside a larger stell ar system, and 
an exponential large scale stellar disc (e.g. Freeman 1970) 
does not conform to modern high accuracy observations, in 
particular to photometric data. With the surface photom- 
etry reaching low surface brightness limits, down to 27- 
28 mag arcsec -2 in the r-band, the brightness profiles of 
most large scale stellar discs cannot be fitted by a single 
component exponential law: 90 pe r cent of discs turn to b e 
either truncated or antitruncated (|Pohlen fc Truiillolliooej) . 
Bulges also appear to be far from de-Vaucouleurs' spheroids: 
when approximated by a Sersic law with a free power param- 
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eter, they show surface brightness p rofile shapes spanning a 
range of n, with the mean n ~ 2 ([And rcdakis et al.l 199a; 
Seigar fc Jamesl Il998l : ICrahaml l200ll : iMollenhoff fc Heidtl 
20011) . Som e morphological types , nam ely, spirals later 



than Sbc ([An drcdakis & Sander j 19941) a nd lenticulars 
ijBalcells et al.ll2003l : lLaurikainen et al.ll2005l h exhibit expo- 
nential surface brightness profiles for their bulges. The scat- 
ter in the concentration parameter n is thought to be re- 
lated to variety of bulge formation mechanisms; a range of 
scenarios restricts also dynamical properties of the bulges of 
different types. 

Currently, all bulges are divided into two categories, 
"classical" high luminosity bulges which are suggested to 
be formed by fast violent events like mergers, and "pseu- 
dobulges". The bulg es of the l atter type promulgated 
by J. Kormendy (see iKormen dv 1993 and earlier confer- 
ence contributions), are thought to form by secular evo- 
lution from the gaseous and stellar material of the disc 
(|Kormendv fc Kennicuttl |2004 ). in particular through the 
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central disc heating by bars |Combes fc Sanders! Il98ll ). 
hence they hav e to resemble di s cs by some of their dynamical 
characteristics. iBalcells et all (|2007t > noted a presence of a 
lot of separat e discy stellar struct ures in the centers of spiral 
galaxies, and lErwin et all ||2003h demonstrated that some- 
times inner discs in early-type disc galaxies might resem- 
ble bulges. Although, following this fashion, a great variety 
of centr al structures often foun d in early-type disc galaxies 
(see e.g. lErwin fc Sparkdliool ) is now trea ted as the pseu- 
dobu lges including nuclear discs and bars (|Drorv fc Fisherl 
2007). we would like to stay on the classical point of view 
that the principal difference between discs and bulges is re- 
lated to their thickness: discs are thin, flat, roughly two- 
dimensional structures, and bulges must be thick and three- 
dimensional by definition. Pseudobulges being produced by 
secular-evolution m echanisms may have e xponential surface 
brightness profiles |Friedli fc Bend 1 19951 ) , but they differ 
principally from the exponential stellar discs by their thick- 
ness. To put a quantitative criterion to distinguish between 
discs and spheroids, let us take a look on the phenomenology. 
E. Hubble classified elliptical galaxies by their shape ranged 
between E0 (the axis ratio of 1) and E7 (the axis ratio of 
3); however many former E7 galaxies are now thought to be 
lenticulars, NGC 3115 being a famous example. As for the 
large-scale stellar discs, the closest axis ratio to the transi- 
tion towards spheroids is perhaps demonstrated by the thick 
disc of the edge-on spiral galaxy NGC 891: its ratio of the ex- 
pon ential scalelength to the exponential scaleheight is about 
3.3 |lbata et al.ll2009l ). Hence, we consider that the ratio of 
scalelength to scaleheight of about 3 is a reasonable frontier 
between spheroids and discs. 

While analysing a photometric structure of a disc galaxy 
inclined to the line of sight, we can trace the ellipticity of 
its isophotes to establish a transition radius where the thick 
central structure, the bulge, delegates a dominance to the 
thin disc: it is the radius where the isophote ellipticity stops 
raising. For a galaxy at a small inclination, it is very diffi- 
cult to distinguish between its exponential pseudobulge and 
a disc without a 3D dynamical model of a galaxy. We at- 
tempted to construct such a model for the galaxy NGC 7217, 
for which there is still no agreement on the structure and 
origin of its inner subsystems. 

NGC 7217 is a giant early-type spiral galaxy seen al- 
most face-on. The distance to the galaxy adopted in our 
paper, 1 8.4 Mpc, was e stimated from the Tully-Fisher re- 
lation bv lRusse if (j2002T l. It corresponds to the spatial scale 
of 0.08 kpc arcsec -1 . NGC 7217 is li s ted in the catalogue 
of isolated gal axies by Karachcnts eval ()l973l ): in a more re- 
cent study bv lBettoni et al. (|200ll ). the density of its envi- 
ronment is also estimated as zero. The galaxy whose evo- 
lution is supposedly free of the environmental influence, 
demonstrates a set of enigmatic structures which can be 
best explained by (a set of) minor mergers. First of all, it 
has three starforming rings , at radii of 11, 33, and 75 arc- 



sec (jButa fc Crocker! Il993l ) looking like resonance struc- 
ture s requiring the presence of a n on-axisymmetric poten- 
tial (jVerde s- Monten egro et aDll995l ). while the galaxy itself 
is unbarred. iButa et al.l (1 19951 ) proposed to decompose the 
whole galaxy into a small disc and a large mildly triaxial de- 
Vaucouleurs spheroid dominating at all distances from the 
center. The dominance of the bulge triaxial potential over 
the whole galaxy caused the formation of three resonance 



rings with high density gas concentration and ongoing star 
formation. Later we undertook our own decomposition of the 
NGC 7217 structure and have found the two large exponen- 
tial components with different scalelength^], 12.5 arcsec, or 
1 kpc for t he inner one and 35.8 arcsec, o r about 3 kpc for the 
outer one (|Sil'chenko fc Afanasievll2000h . The inner compo- 
nent has larger intrinsic ellipticity than the outer one and so 
may be oval providing the necessary triaxiality of the poten- 
tial to put the rings at the resonance radii. The bulge, if any, 
is small, confined within the innermost (nuclear) ring, and 
exponential. To choose between the alternate ways of decom- 
position, the kinematical data and a dynamical model are 
needed. The situation is complicated by the fresh (?) rem- 
nants of a minor merger presented by counte rrotating stars 
in the inner disc (|Merrifield fc Kuiikenll 19941) and the inner 
gas p olar disc within R ~ 4 arcsec (ISil'chenko fc Afanasievl 

In this paper we present the results of long slit spec- 
troscopy of NGC 7217. Our goal is the diagnostics of the 
nature of the two large scale exponential structures in this 
galaxy. The inner exponential structure may be an inner 
disc (a thin structure) or a pseudobulge (a thick structure). 
Since NGC 7217 is seen nearly face-on, the isophote elliptic- 
ity is a poor indicator of the stellar component thickness as 
it is close to zero over the whole galaxy. But with the stel- 
lar velocity dispersion profiles along the major and minor 
axes we can try to estimate the scaleheights of the stellar 
components from dynamical considerations. In Section 2 we 
describe our observations, data reduction and analysis, in 
Section 3 we provide our new estimates of the stellar popu- 
lation properties at different radii. The principal Section 4 
gives our dynamical consideration of the NGC 7217 large 
scale structures, and Section 5 contains a discussion about 
the origin of the galaxy structure including the compari- 
son with numerical simulations. In Section 6 we give a brief 
summary of our results. 



2 OBSERVATIONS AND DATA REDUCTION. 

The spectroscopic observations were carried out with the 
SCORPIC0 universal spectrograph (jAfanasiev fc Moiseevl 
l2005h installed at the prime focus of the Russian 6-m 
Bol'shoy Teleskop Azimutal'nyy (BTA) operated by the Spe- 
cial Astrophysical Observatory, Russian Academy of Sci- 
ences. We used the VPHG2300G grating providing an in- 
termediate spectral resolution (R ~ 2200) in a relatively 
narrow wavelength region (4800 < A < 5500 A) however 
containing a rich set of strong absorption line features mak- 
ing it suitable for studying both internal kinematics and stel- 
lar populations of a galaxy. The chosen spectral range also 
includes several emission lines, H/3, [Om], and [Nl], which 
we used to derive the gas kinematics and line ratios. The 
slit was 1.0 arcsec wide and 6 arcmin long. The 2kx2k EEV 
CCD42-40 detector used in the 1x2 binning mode provided 
a spectral sampling of 0.37 A pix -1 and a spatial scale of 
0.357 arcsec pix -1 . 



1 The bulge also d emonstrates an exponential pr ofile with a scale- 
length ~ 4 arcsec llSil'chenko fc Afa nasicv 200^). 

2 For a description of the SCORPIO instrument, see 
http : //www . sao . ru/hq/moisav/scorpio/scorpio .html 
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We observed NGC 7217 in two slit positions going 
through the centre, along minor and major axes of its in- 
ner isophotes. The major axis spectrum, P.A.=81 deg, was 
obtained on 6/Oct/2008 with an integration time of 80 min 
under atmosphere conditions with good transparency and 
intermediate image quality of 1.7" FWHM corresponding to 
the median seeing at the telescope site. The minor axis data 
at P.A.=169 deg were collected during two nights October, 
6, 2008 and October, 8, 2008 with a total exposure time of 
75 min under bad variable transparency and seeing qual- 
ity of ~ 3 arcsec making them notably shallower than the 
major axis spectra. We obtained the following calibrations: 
night time internal flat field and arc line spectra, the GD 248 
spectrophotometric standard star, and high signal-to-noise 
twilight spectral frames. 

Data reduction and analysis for the spectral data of 
NGC 7217 was identic al to that of the l enticu lar galaxy 
NGC 6340 presented in IChilingarian et all |2009l ). We refer 
to that paper for all details, here we give only essential in- 
formation required for understanding our data analysis and 
interpretation. 

The primary data reduction steps included bias sub- 
traction, fiat fielding and cos mic ray hit remova l using the 
Laplacian filtering technique l|van DokkumllioOlh . Then, we 
built the wavelength solution by identifying arc lines and fit- 
ting their positions using a 3rd order two-dimensional poly- 
nomial along and across dispersion, and linearised the spec- 
tra. The obtained wavelength solution had fitting residuals 
of about 0.08 A RMS. 

The SCORPIO spectrograph has significant varia- 
tions of the spectral line-spread-functi on (LSF) along and 
across the wavelength direction (see iMoiseevl 120081 and 
IChilingarian et al ] |2009l for details). In our observations of 
NGC 7217, we used the peripheral regions of the slit beyond 
2 arcmin from the galaxy centre to estimate the night sky 
spectrum, therefore it was very important to take the LSF 
variations into account in order to minimize the effects of the 
sky subtraction artefacts on the data analysis. We mapped 
the LSF by fitting the high-resolution (R = 10000) Solar 
spectrum against the twilight spectra at 64 positions along 
the slit in five slightly overlapping wavelength segments cov- 
ering the spectral range of the SC ORPIO setup with the 
pena lized pixel fitting technique (|Cappellari fc Emselleml 
120041 ). We us ed the Gauss-Hermite param etrization up to 
the 4th order llvan der Marel fc Franxlll993h to represent the 
LSF shape. 

Then, we modelled the night sky spectrum at every po- 
sition al ong the slit follow i ng th e procedure described in 
detail in IChilingarian et alJ (|2009h where the main step was 
the parametric signal recovery applied to the sky spectrum 
taken in the outer regions of the slit in the Fourier space as 
follows: 



f(x,X) = F-\F(f (sky, A)) 



F(£(sky)) 



(1) 



where f(x, A) denotes a sky spectrum at the position x along 
the slit with its parametrized LSF C(x); /(sky, A) is the 
night sky spectrum in any region of the slit with the LSF 
£(sky), and F, F^ 1 are the direct and inverse Fourier trans- 
forms respectively. The night sky model created in this fash- 
ion results in a nearly Poisson quality of the sky subtraction 



which is absolutely crucial for our analysis of the external 
regions of the galaxy. 



3 STELLAR POPULATIONS AND INTERNAL 
KINEMATICS OF NGC 7217 

We derived the parameters of internal kinematics and stel- 
lar popul ations of NGC 7217 by fitting high-resolution pe- 
GASE.hr (|Le Borgne et all [20041 1 simple stellar population 
(SSP) models against ou r spectra with the NBur sts full 
spectral fitting technique (|Chilingarian et al.ll2007bllal) . The 
SSP models were computed using the ISalpete rl d 19551 ) stel- 
lar initial mass function. The fitting algorithm is picking 
up a template from a grid of stellar population models in 
the age-metallicity space convolved with the instrumental 
response of the spectrograph (see previous section), broad- 
ening it with the line-of-sight velocity distribution (LOSVD) 
of a galaxy represented by the Gauss-Hermite parametriza- 
tion up to the 4th order, i.e. v, a, /13, and /14. The models 
are multiplied pixel by pixel by the n order Legendre poly- 
nomial continuum which is used to account for possible im- 
perfections of the flux calibration and for the internal dust 
extinction in a galaxy. All kinematical and stellar population 
parameters are determined in a single non-linear minimiza- 
tion loop therefore reducing possible degeneracies. The \ 2 is 
penalized towards purely Gau ssian LOSVD as explained in 
ICappellari fc Emselleml (|2004l ') in order to stabilise the solu- 
tion in case of low signal-to- noise ratios and/or insufficient 
spectral sampling of the data. We stress that for the dy- 
namical analysis presented below, we use the pure Gaussian 
LOSVD parametrization without x 2 penalization. 

The pegase.hr models are constructed from empirical 
stellar spectra in the Solar neighbourhood where the value of 
a-element abundances is known to correlate with the over- 
all metallicity. Therefore, at intermediate and high metal- 
licities, our SSP models ar e representative of t he so lar a/Fe 
abundance ratio. However. IChilingarian et alJ (|2008l ) showed 
that supersolar a/Fe ratios bias neither age nor metallic- 
ity estimat es when using the NBursts sp ectral fitting tech- 
nique, and lSiPchenko fcA fanasicv (200^) found nearly solar 
a/Fe in NGC 7217. 

NGC 7217 possesses quite strong emission lines in its 
spectrum especially in the regions corresponding to the rings 
with ongoing star formation. Therefore we had to exclude 
from the fitting procedure several narrow 20 A-wide re- 
gions around emission lines (H/3 A = 4861 A, [Om] A = 
4959, 5007 A, and [Nl] A = 5197, 5201 A) reds hifted accord- 
ing to the line-of-sight velocity of NGC 7217. IChilingarian! 
|2009T ) demonstrated that H/3 contains 20 per cent of the age- 
sensitive information at maximum when using the NBursts 
technique in a spectral range similar to ours, therefore ex- 
cluding it from t he fit neither biases age e stimates (see also 
A ppendix A2 in IChilingarian et al. I l2007bl and Appendix B 
in IChilingarian et al" 20081 ). nor degrades significantly the 
quality of the age determination. 

We fitted Gaussians pre-convolved with the SCORPIO 
LSF into the [Om] and H/3 emission lines in the residuals 
of the stellar population fitting and determined the ionised 
gas kinematics independently in the two lines as well as the 
emission line flux ratios. 

The obtained radial distributions of SSP-equivalent 
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Figure 1. Radial distributions of the mean age (top) and mctal- 
licity (bottom) of the stellar populations measured along the ma- 
jor (PA = 81°) and minor (PA = 169°) axes. The vertical dashed 
lines indicate the conventional boundaries of the bulge and the 
inner exponential component set at four scalclcnghts of each com- 
ponent. 



age and metalliticity o f NGC 7217 are shown in Fig. [TJ 
ISil'chenko fc Afanasievl (|2000h analysed the large scale 
structure of NGC 7217 and identified three galaxy compo- 
nents, a bulge dominating at radii 5 — 15 arcsec (0.4-1.2 kpc), 
and two exponential structures interpreted as discs. The in- 
ner and outer discs dominate in the galaxy light profile at 
20 - 50 arcsec (1.6-4 kpc) and 60 - 110 arcsec (4.8-8.8 kpc) 
respectively. 

The age profile of NGC 7217 presented in Fig. Q] ex- 
hibits specific features at the radii dominated by these three 
substructures. In the bulge dominated region located mostly 
inside the nuclear starforming ring (R = 10 — 12 arcsec), the 
mean stellar age decreases from 10-13 Gyr in the centre to 
5 Gyr outwards. In the inner exponential component the age 
stays nearly constant at about 5 Gyr; individual estimates 
at the radii between 20 and 50 arcsec have rms of 0.7 Gyr 
in the major-axis profile and 0.9 Gyr in the minor-axis one. 
The outer exponential component (R — 60—110 arcsec) 
has an intermediate age of about 2-3 Gyr, whereas in the 
broad outer starforming ring (R — 70 — 80 arcsec) it drops 
down to 1 Gyr. Surprisingly, at these radii we also see a drop 
in the SSP-equivalent stellar metallicity down to —0.65 dex 
prominent in the eastern side of the galaxy but hardly de- 
tected in the western side, which is probably connected to 
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Figure 2. Radial distribution of the line-of-sight velocities of 
stars and ionized gas along the major (PA = 81°) axis. The [Oio] 
data are presented only for \R\ < 10 arcsec. The vertical dashed 
lines indicate the conventional boundaries of the bulge and the 
inner exponential component (see Fig. [1} • 



the morphology and location of individual starforming re- 
gions crossed by the slit. The mean metallicity of the bulge 
is close to the solar value, [Fe/H]=— 0.06 dex, while the outer 
exponential component ([Fe/H] =— 0.45 dex) is rather metal 
poor for such a massive galaxy. The inner exponential com- 
ponent has a mean metallicity of about —0.2 dex with a 
notable metallicity gradient around —0.3 dex per dex in ra- 
dius. 

Hence, the existence of the three component s found in 
the photometric data by ISil'chenko fc Afanasievl (|2000l ). a 
bulge and two [probably] exponential discs, is supported by 
specific features in the stellar population profiles. Later, we 
will propose the evolutionary scenario of NGC 7217, where 
the two discs may have different origin. 

We present the derived major axis (P.A.=81 deg) pro- 
files of line-of-sight velocities for the gaseous and stellar 
components in Fig. [2] The H/3 emission line component re- 
mains clearly visible along the whole extent of the galaxy 
where we can measure the stellar kinematics. The [Om] is 
much weaker than H/3 at radii R > 10 arcsec. The veloc- 
ities obtained from the two emission lines agree well, but 
those determined from H/3 in the outer region have 4-5 
times smaller uncertainties. Therefore, we constructed the 
combined gaseous line-of-sight velocity profile displayed in 
Fig. [2] from [Oin] and H/3 kinematics at radii below and 
above 10 arcsec correspondingly. Both, stellar and gaseous 
velocity profiles have similar shapes including even some 
particular features such as a notable drop at R ~ 50 arcsec. 
They are also very symmetric. Therefore, for the dynami- 
cal analysis presented in the next section, we folded them 
along the galaxy centre and averaged the values of rotation 
velocities from both sides. The gaseous rotation velocities 
(~ 150 km s _1 ) exceed the stellar ones (~ 120 km s" 1 ) 
by nearly 30 km s _1 at all radii, which is an observational 
manifestation of the asymmetric drift. The difference stays 
nearly constant between 10 and 70 arcsec suggesting similar 
dynamical statuses of the inner and outer exponential stellar 
structures. These structures are probably discs, but quite 
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Figure 3. Radial distributions of the line-of-sight stellar velocity 
dispersion measured along the major, PA = 81°, (top) and minor, 
PA = 169°, (bottom) axes. The vertical dashed lines indicate the 
conventional boundaries of the bulge and the inner exponential 
component (see Fig.[TJ. 



high asymmetric drift implies that they should be rather 
thick. 

In Fig. [3] we show the observed radial distributions of 
the line-of-sight stellar velocity dispersion along minor and 
major axes of NGC 7217. Since the galaxy is seen nearly 
face-on, the vertical (a z ) component of the disc velocity el- 
lipsoid creates the main contribution to the observed veloc- 
ity dispersion. However, there is a notable difference between 
major- and minor-axis velocity dispersion profiles suggest- 
ing significant contributions of the tangential (cr v ) and ra- 
dial (<jr) velocity dispersions to the values along the major 
and minor axes respectively. The main feature of the stel- 
lar line-of-sight velocity dispersion profiles of NGC 7217 is 
a minimum in the area of the inner exponential component 
and a smooth raise outwards. This behaviour is quite unex- 
pected: as we demonstrated, the outer stellar component is 
younger than the inner one, so it would be more natural if 
the stellar subsystem formed recently from the dynamically 
cold gas is itself dynamically cold. However, the velocity dis- 
persion rise beyond R « 50 arcsec is statistically significant 
and seen in both minor- and major-axis velocity dispersion 
profiles. In Section 5 we will give the possible explanation 
of these features. 



4 DYNAMICAL ANALYSIS: TWO VERY 
DIFFERENT DISCS IN NGC 7217 

We used our kinematical data to recover the stellar velocity 
ellipsoid and reconstruct the radial distributions of all stel- 
lar velocity dispersion components. Then, we calculated the 
stellar disc thickness profile from the a z radial distribution. 
By using relations describing the disc equilibrium, we did it 
independently from the major- and minor-axis kinematical 
data. 

For an intermediately inclined galaxy, the two in-plane 
components of the stellar velocity dispersion, or and cr v , are 
related to the observed minor- and major-axis line-of-sight 
velocity dispersions according to the following relations: 



CTfos,mm(-RcOS l) 



2 ■ 2 • i 2 2 • 

o R sin i + <7, cos i , 



los.maj 



(R) 



(2) 



They also include the contribution from the vertical com- 
ponent of the velocity ellipsoid o z . Thus, by measuring the 
variations in the line-of-sight dispersions along the princi- 
pal axes, we can obtain only a linear combination of <tr{R), 
o lp (R) and o z (R). Therefore, we need additional information 
in order to derive all three velocity dispersion components 
from these equations. To close the system of equations ([2]), 
we can use some dynamical relations, which are valid if the 
system is in equilibrium. One such relation connects the ve- 
locity dispersion component s or(R) and o v (R) via the mean 
azimuthal velocity of stars (|Binnev fc Tremainelll987t ) 



4 = i|i + 



dlnv v 
dlnR 



(3) 



If the rotational velocities are unknown, we can use the local 
circular speed of gas v c instead of v v in order to relate o R 
and at 



5( 1 + 



<91n« c 
dlnR 



(4) 



This relation is true if most orbits in a disc are quasi-circular. 

Adding any of these relations, Eq. [3] or Eq.[4] to the two 
presented above, enables us to recover radial distributions 
of all three velocity dispersion components. For the first 
time this technique w as applied to the data for NGC 488 by 
iGerssen et al.l (1 19971 ). However, this procedure creates very 
noisy output when applied directly to the data, because it 



los.maj 



with very close values, as well as the numerical 



the kinemati cal profiles and to find the best fitting solu - 
tion (see e.g. IGerssen et alii 19971 , l200d : IShapiro et al]|2003l ) 
which will, however, depend on the adopted parametriza- 
tion. For thi s reason, we adopted the less parametric ap- 
proach as in iNoordermeer et al.l (|2008h . We approximated 
all kinematical profiles using polynomials and calculated all 
quantities including their derivatives analytically. But even 
then, the subtraction of the major axis velocity dispersion 
profile from the minor axis one results in unreliable and 
ambigious solution. To avoid this, we use the asymmetric 
drift equation. The deprojected major-axis gas rotation is 
a measure of the circular speed v c , while the deprojected 
major-axis stellar velocities allow us to determine the mean 
rotational motion v v related to the radial velocity dispersion 
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or. Therefore, we can obtain the or profile using the major- 
axis velocity p rofiles for gas and stars an d the asymmetric 
drift equation dBinnev fc Tremaindlf 98?l ): 



-2 2 
V v = CTr 



- 1 - 



<91nS _ ainqj 
dlnR ~ dlnR 



(5) 



where E is the stellar surface density. Provided that the 
stellar mass-to-light ratio stays nearly constant ((M/L)i — 
1.86 in the Solar units from the PEGASE.2 models, see 
iFioc fc Rocca-Volmerangel 1 19971 ). we can use the surface 
brightness instead of surface density E in Eq. [5] We assume 
that the /-band photometric data trace old stellar popu- 
lation whose internal kinematics we are studying. Because 
of the logarithmic derivative, the exact choice of the M/L 
ratio is not critical. The asymmetric drift equation can be 
used directly to obtain the radial velocity dispersion profile, 
provided that the ratio between tr^/o- R is determined from 
Eq. [3] or Eq. [4] However, Eq. [5] requires the knowledge of 
the or radial gradient. We assume that in the outer regions 
of the galaxy between 30 and 70 arcsec, or may decline ex- 
ponentially as tx exp(— R/hkin) and use the last term in the 
form 2i?//ikin, where /ikin is a free parameter. 

The choice of a radial velocity dispersion profile expo- 
nentially declining with the radius is physically motivated. 
There is a conventional assumption that the disc thickness 
and the z-component of the velocity dispersion a z are con- 
nected via the vertical e quilibrium condition for an isother- 
mal layer (|Spitzerlll942h 



a z (R) =ttGZ(R) Zq, 



(6) 



where zo is the half-thickness of a homogeneous layer. For a 
mass-to-light ratio constant with radius, constant zq and an 
exponential brightness profile it yields: 



cr z (R) oc exp(— R/h) . 



(7) 



where h is a disc scalelength. 

The disc heating theory suggests that the range 
in the velocity anisotr opy o z /<jr, is about 0.4 — 0.8 
d Jenkins fc Binnevlll990l). which is consistent with observa- 



tions. 

modelled the data for several galaxies and estimated the 
values for the velocity anisotropy a z /<7 R between 0.5 and 
0.7, s imilar to the solar neighbourhood (|Dehnen fc Binnevl 
1 19981 ). 

Since the ratio (j z /ctr is close to constant, Eq. [7]yields: 



ctr oc exp(— R/2 h) . 



(8) 



It implies th at fekin = 2 h. This app roach does work for 
our Galaxy (| Lewis fc Freeman] Il989l ) but in the case of 
NGC 7217 we did not fix the value of /ikm and estimated 
it iteratively by fitting the radial velocity dispersion profile 
with the exponential law and substituting the value of ftkin 
in the asymmetric drift equation. We needed this value to 
estimate only the last term in Eq. [5] The exact value of 
this term may affect the final results but the effect it quite 
insignificant. We were changing the range of R used to fit 
or by exponential law and did not notice any difference be- 
tween final results. Having derived the or radial distribution 
from the Eq. [5] we can then compute the radial profile of 
the vertical component of the velocity dispersion a z using 
one of the Eq. [2] We choose to use both of them to control 
the reliability of our results. 
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Figure 4. Radial distributions of the radial and vertical velocity 
dispersions reconstructed from the major- and minor-axis kine- 
matical profiles. Full lines - profiles calculated via the Eq.|3] sym- 
bols — profiles obtained by using the Eq. [4] Error bars are shown 
only for profiles indicated by solid lines. The vertical dashed lines 
indicate the conventional boundaries of the bulge and the inner 
exponential component (see Fig.[TJ. 



Finally, we can obtain the thickness profile from Eq. [6] 
which is valid for an isothermal layer. The assumption about 
constant velocity dispersion along the z-direction for discs 
is in agreeme nt with the results of iV-body simulations (see 
e.g. Fig. 2 in lSotnikova fc Rodionovll2006l ). But one should 
keep in mind that Eq. [6] gives the upper limit for the disc 
thickness. If there is a massive dark halo, the disc thickness 
will be lower for the same value of a z . 

In order to deproject the velocity profiles we need 
to know the inclination of the rotation plane to the line 
of sight. NGC 7217 has a roundish appearance, how- 
ever quite high observed (i.e. projected) rotation veloc- 
ities suggest that is not seen completely face-on. The 
isophote ellipticity analysis assuming the infinitely thin 
stellar disc provided the dis c inclination estimates in the 
28 deg ( Noordcr meer fc van der Hulst 



2007 ; IVerdes-Montenegro et al.lll995l ; ISanchez-Portal et al 
20001 ); whereas the kinematical analysis usuall y results in 
slightly higher values i = 3 — 31 deg (see iButa et al.l 



19951; I Noordermeer et ah! 12 005 for the Hi map analysis and 
Sil'chenko fc Moiseevll2006l for the stellar kinematics) . Keep- 



ing in mind possible considerable thickness of the stellar disc, 
we adopt the inclination value of i — 30 deg. Varying i be- 
tween 26 and 35 deg changes the resulting disc thickness 
estimates by no more than 10 per cent. 

We present the reconstructed radial distributions of ctr 
and a z in Fig. [4] We estimated unc ertainties using a boot- 
strapping method |Press et al.ll 19921 ) . We computed or and 
o"z assuming that the errors of observational quantities used 
in the procedure are distributed according to the normal law 
and changing input data points according to these errors. We 
made several thousands realisation of the procedure, aver- 
aged resulting values, and estimated the dispersion of all the 
simulations at each point. 

The data in the Fig. [4] suggest that the dynamically 
coolest subsystem of NGC 7217 as concerning the verti- 
cal velocity dispersion is the inner stellar exponential com- 
ponent. Therefore it is certainly a disc and not a pseudo- 
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Figure 5. Relation between the ratios of observed maximal ro- 
tation velocities and velocity dispersions and isophote ellipticities 
for early-type galaxies, bulges (open symbols) and pseudobulgcs 
(filled symbols). The models of prolate and oblate isotropic ro- 
tators are overplotted by dashed lines. NGC 7217 is shown with 
a red cross being clearl y abov e the locus of pseudobulges from 
iKormendv fc Kennicuttl l |2004h . 



bulge. Indeed, in Fig. [5] we plot the inner exponential com- 
ponent of NGC 7217 at the diagram of Binney-Kormendy. 
With its e ~ 0.1 and Dmax/fios = 1-24 it stays well above 
all the pseudobulges dynamically confirmed by Kormendy. 
At R > 30 arcsec (>2.4 kpc) the vertical velocity dis- 
persion component slowly raises along outwards reaching 
110 km s _1 at the location of the outer starforming ring. 
The an profile decreases along the radius as expected, but 
in the transition region where the outer disc starts to dom- 
inate (R — 40 — 50 arcsec; 3.2-4.0 kpc), we see a break 
of the monotonic fall. The whole an profile can be divided 
into two parts corresponding to the two exponential stellar 
components. 

The disc "thickness profiles" at radii R = 15 ... 80 arc- 
sec reconstructed from the major- and minor-axis kinemat- 
ics are shown in Fig. [6] While the qualitative agreement 
between them exists, quantitatively they diverge at some 
radii. As it was mentioned in Section 2, the minor axis data 
are of a worse quality than the major axis ones but they 
are featureless as concerning the line-of-sight stellar veloc- 
ity dispersion. In this sense, both datasets are worth each 
other. One of the sources of the profile discrepancy might be 
the adopted inclination angle and the discrepancy between 
measured distances along the major axis and deprojected 
distances along the minor axis. Similarly to the a z profile 
behaviour, we see that the inner stellar disc is relatively thin 
(zo = 0.2 . . . 0.7 kpc). Hence, at its inner boundary the ratio 
of its scalelength to the scaleheight is about 4-5. The outer 
disc flares vigorously reaching the half-thickness of 2.5-3 kpc 
at the outer starforming ring position. 
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Figure 6. Radial distribution of the stellar disc scaleheight - the 
reconstructed "edge-on" view. Solid lines are the profiles calcu- 
lated via Eq. [3] while symbols are for the profiles obtained by 
using Eq. [4] Error bars are shown only for profiles indicated by 
solid lines. The vertical dashed lines indicate the conventional 
boundaries of the bulge and the inner exponential component 
(see Fig. Q}. 



5 DISCUSSION 

5.1 Explaining the structure of NGC 7217 

From the analysis of the stellar and gaseous kinematics (line- 
of-sight velocities and velocity dispersion radial distribu- 
tions) over an extended area of the giant early-type disc 
galaxy NGC 7217, we reconstructed the velocity ellipsoid 
and calculated the stellar disc scale heights at radii from 20 
to 70 arcsec (1.6 to 6 kpc). Earlier. ISil'chenko fc Afanasievl 
120001) demonstrated that this radial range contains two stel- 
lar substructures having exponential surface brightness pro- 
files with different scalelengths, the inner one with relatively 
low h, and the outer one with more extended (standard for 
giant galaxies) scalelength. Our analysis reveals that the in- 
ner component has the scaleheight of 0.2-0.7 kpc and the 
outer one about 1-3 kpc. Hence, we can classify the inner 
and outer components as thin and thick stellar discs respec- 
tively. 

Having applied the full spectral fitting NBuRSTS tech- 
nique with a grid of high-resolution SSP models, we de- 
duced mean ages and metallicities of stars. The inner disc 
exhibits the intermediate age of about 5 Gyr and the mean 
metallicity of —0.2 dex with a strong negative metallicity 
gradient along the radius. The outer disc is quite metal- 
poor ([Fe/H]= —0.4 dex) and relatively young (2 Gyr). The 
outer disc harbours a prominent starforming ring with the 
mean age of about 1 Gyr at the radius of R ~ 75" very 
well visible in UV data from the GALEX satellite and 8/im 
NIR data from the Spitzer Space Telescope. Interestingly, 
in this ring the mean stellar metallicity falls down to a very 
low value of —0.7 dex. The galaxy nucleus and the bulge of 
NGC 7217 inside the nuclear starforming ring are very old 
(~ 13-15 Gyr). 

How can we explain the complex structure of 
NGC 7217? In principle, the fact that the outer stellar disc 
is younger than the inner one is in line with the paradigm 
of the 'inside-out' disc formation. But the sharp boundary 
between the two discs, where the properties including stellar 
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Figure 7. Ratio between observed gas density and critical den- 
sity according to iKenni cutj (fl 989l^ and the two-fluid instability 
llJog fc Solomonlll984lEfstathioull2000h criteria. The dashed line 
indicates the threshold for star formation. The H i surface density 
profile was taken from INoordermeer et al, ] ||2005| 1. 



age and metallicity as well as the radial and vertical scale- 
lengths of the star distribution change abruptly, hints to a 
temporal gap between the inner and outer disc formation, 
or to some catastrophic event that provoked the expansion 
of star formation into the outer ring area. 

Interesti ngly, the mol e cular gas concentrates in the 
inner disc (|Combes et all |2004 ) while the neutral hy- 
drogen is__ol2served_^n__the_narrow r ing in the outer 
disc dVerdes-Montenegro et all 1 19951 ; iButa et al] 1 19951 ; 
iNoordermeer et al .1 120051 ). However, star formation is almost 
absent in the inner disc (though it burns in the nuclear 
starforming ring between the bulge and the inner disc) an d 
is very prominent in the outer disc (|Battinelli et al.ll2000l ). 
INoordermeer et al] (|2005l ) estimated the gravitational sta- 
bility of the gas in the outer ring a nd the thre s hold d ensity 
for star formation according to the iKennicu tt (1989) crite- 



(9) 



3.36G ' 

where k is the epicyclic frequency, that can be derived from 
the rotation and c is the gas velocity dispersion, which can 
be assumed as a constant value of 6 km s _1 . For the dimen- 
sionless quantity a, Kennicutt derived a value of 0.67 from 
the empirical study of star formation cutoffs in spiral galax- 
ies. Dynamically, it means that one should take into account 
non-axisymmetric modes while considering t he stability of a 
disc, because for ax i symm etric modes a = 1 |Toomrelll96"3 ). 
INoordermeer et al] (|2005l ) demonstrated that even in the 
middle of the Hi ring of NGC 7217, the gas density does 
not exceed 25 per cent of the critical one. The estimate by 
INoordermeer et al] (2005) did not take into account the exis- 
tence of a stellar disc. Having restored the velocity dispersion 
profile or for stars, we can evaluate the gase ous disc stability 
in the framework of th e two-fluid approach (|jog fc Solomon! 
11984 lEfstathioull2000l ) 
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Figure 8. Minor-axis kinematics and emission line properties of 
NGC 7217 from the analysis of the H/3 (red, diamonds), [Olll] 
(A = 5007 A, blue, triangles), and [Nl] (A = 5197/5200 A, green, 
squares) lines in its spectra with the best-fitting stellar population 
models subtracted. The panels show (top to bottom): radial ve- 
locities, velocity dispersion, line fluxes, logarithms of the emission 
line ratios: [Om]/H/3 (blue) and [Nl]/H/3 (green). 



where a and b are the ratios of stellar to gas velocity dis- 
persion and surface densities respectively, and g(a, b) is a 



function derived numerically bv lEfstathioul lj2000h . This ap- 
proach results in the critical den sity estimate at least twic e 
lower than the value obtained by INoordermeer et al] (|2005h . 
But even now the critical density remains too high (see 
Fig. [7}. Therefore, the gas in the outer ring must be gravi- 
tationally stable, and why the star formation proceeds there 
is still an open question. The very existence of a system of 
starfor ming, gas-populated rings in an unbarred galaxy is a 
puzzle. ICombes et al] (|2004l ) discussed a possibility of a past 
bar, strong or weak, which had formed the rings and then 
had dissolved after the gas inflow had supplied mass into 
the center. Now we see that this scenario experiences strong 
difficulties: the stellar population in the NGC 7217 centre 
is very old, and it is clear that no significant star formation 
took place there for the last 5 Gyr in order to provide mass 
concentration and the subsequent bar destruction. 

5.2 Inner polar disc 

In Fig. [8] we present the minor-axis kinematics and emis- 
sion line properties in the central region of NGC 7217 de- 
rived from the analysis of the stellar population spectral 
fitting residuals. In the nuclear region the [Nl] doublet is 
clearly detected in addition to the oxygen and hydrogen 
lines. The [Olll] kinematics is evident of either a rotating 
disc or an outflow from the galaxy centre which can be in- 
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Figure 9. Inner polar disc in NGC 7217 as seen with the HST. 
This false colour composite is made of FA50W, F606W, and 
F814W WFPC2 images with the Ha fluxes from the ACS F658N 
images added to the red channel using the iLupton et al, | 1120041 1 
algorithm. The strongly inclined gaseous structure seen here in 
the centre is a polar disc because the outer disc is seen almost 
face-on with its major axis nearly horizonthal in this plot. 



duced by an active nucleus. We measured the emission line 
ratios (see bottom pan el in Fig. [8} and used the itera tool 
(I Groves fc Alien! |201Ch providing a large collection of emis- 
sion line models for different excitation mechanisms to per- 
form the diagnostics. The extremely high log([Ni]/H/3) > 
ratio given moderate log([Om]/H/3) ~ 0.45 is consistent 
with either a shock ionisation for solar or subsolar metal- 
licities of the ISM, or an AGN having a significantly super- 
solar metallicity (e.g. > +0.3 dex). Given very moderate 
metallicity of the stellar population in the central region of 
NGC 7217, the latter possibility can be ruled out at a high 
level of confidence. Therefore, we are probably observing an 
inner polar disc in NGC 7217, which is supported by the 
direct inspection of the HST imaging data in the narrow 
Hq+ [Nil] A6583-ce ntered filter (F i g-Ell. 

According to iKarachentseval |l973T l. NGC 7217 is an 
isolated galaxy, not a member of any g roup, and its environ- 
ment is emp ty of satellites. Howeve r, I Merrifield fc KuiikerJ 
(|l994h and ISil'chenko fc Afanasievl (|2000h suggested the 
presence of a counterrotating stellar component in the cen- 
tral region of NGC 7217; moreover, we see also the inner 
polar gaseous disc. Hence, accretion or minor merger(s) in 
the past seem inevitable. If we assume that the galaxy cen- 
tre is its oldest part then we can suspect the vertical im- 
pact by a gas-rich satellite. Then the 0.3 kpc nuclear polar 
disc in NGC 7217 is a remnant of the satellite gaseous com- 
ponent which conserved its initial orbital momentum. The 
close analogue of such an event is the Sagittarius dSph dis- 
ruption by the Galaxy accompanied by the re-distribution 
of the satellite material along the polar large circle. Polar 
orbits are known to be stable, hence this gas cannot accrete 



to the very centre explaining why the fuel for the star for- 
mation in the centre of NGC 7217 has been absent until now 
evidenced by the old stellar population there. 

Such a vertical impact is also able to produce a system 
of starforming/stellar rings by exciting a running compres- 
sion wave in a large scale galaxy disc (|Athanassoula et alj 
Il997l). Moreover . the si mulations of a Cartwheel-like galaxy 
by Mapelli et al.l (|2008l ) for 1.5 Gyr after the impact clearly 
show that sharp stellar rings formed initially expand over 
a large area during 1 Gyr and form a structure resembling 
a low surface brightness outer stellar disc. The two-tiered 
exponential stellar disc of NGC 7217 could have been built 
by such a minor merger. But we still need a triaxiality of 
the potential in order to produce a counterrotating compo- 
nent in the inner disc during the same event: the polar gas 
in the tumbling triaxial potential has to warp in the outer 
part in such a sense that it ar rives to the rotational plane 
of stars with the opposite spin (jvan Albada et alll982l ). Be- 
ing compressed, this gas would be consumed by star forma- 
tion leaving a counterrotating stellar component. Besides, 
the triaxial potential also creates long-living resonance rings 
in contrast to short-living collisional rings. From our data we 
cannot distinguish between a triaxial dark-matter halo and 
a low-contrast e xtended stellar spheroid (see discussion in 
iButa et al.lll995l ). 



5.3 Insights from simulations 

Simulations of vertical impacts producing ring structures re- 
veal that the vertical oscillations of disc stars excited by the 
intruder' s passage result in stron g thickening of the stel- 
lar disc (|Hernquist fc Weill ll993T ) which we probably ob- 
serve in the outer disc of NGC 7217. The small thickness 
of the inner disc can be perhaps explained by its high den- 
sity increased after the impact by the counterrotating com- 
ponent and also by the presence of the significant molec- 
ular gas content (|Combes et al. I l2004h . Interestingly, in the 
recognized vertical collision product, the famous Cartwheel 
galaxy, the molecular gas concentrates in the central region 
while the neutral atomic hydrogen is confined to the outer 
ring (|Horellou et al.lll998T l being very similar to NGC 7217. 

One can suggest another explanation of the inner 
disc structure. It could have been formed from the ma- 
terial of a small satellit e during its complete disruption 
l|Eliche-Moral et al.ll2006l) . This new inner disc can be older 
than the outer one because the infalling stars were initially 
older. It can also be dynamically cold because its material 
had experienced the orbital circularisation. However, this 
scenario requires a nearly in-plane encounter. The unusu- 
ally high thickness of the outer component and the presence 
of the inner polar disc in NGC 7217 make us to prefer en- 
counters at higher inclinations. 

The minor merger scenario with a vertical impact re- 
quires this event to have occurred long time ago. The galaxy 
needs time to convert the counterrotating gas into stars and 
to establish the mean stellar age of 5 Gyr in the compact 
high-density inner stellar disc. Then the current star forma- 
tion in the outer disc involving very low-metallicity gas does 
not relate to this event and to its collision rings. It may 
be provoked by later gas infall in a tidal tail developped 
during the minor merger. This mechanism can explain the 
patchy and filamentary edges of the Hi ring in the map of 
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NGC 7217 (|Noordermeer et alj [20051 ) and the presence of 
the shock-induced star formation in the low-density gas ring 
which had to be gravitationally stable under the quiescent 
conditions. 

We explored intermediate resolution (0.2 kpc) TreeSPH 
simulations of minor mergers prov ided by the GalMer 
databas^fl (|Chilingarian et all |2010| ). Presently available 
simulations include interactions of a gas free giant lentic- 
ular galaxy ("gSO") with a 10 times less massive dwarf 
galaxy. The morphology of a dwarf spans the entire Hubble 
sequence from non-rotating ellipticals ("dEO") to bulgeless 
discs ("dSd"). Here we consider only gSO-dSd mergers, as 
this configuration has the largest quantity of gas compared 
to other morphologies of a satellite galaxy. The inclination of 
the orbital plane to the rotation planes of the two interacting 
galaxies is 33 deg for a giant SO and 130 deg for dSd. The 
two sets of numerical experiments include interactions on 
prograde and retrograde orbits with different initial orbital 
momenta and motion energies (24 in total). The simulations 
were run for a total duration of 3 Gyr. 

The bar is always formed in the gSO stellar disc. Merger 
remnant morphologies vary quite a lot, but all of them can 
be classified into two families related to the mutual orien- 
tation of the infalling satellite internal angular momentum 
with respect to that of the orbital motion. 

All retrograde encounters result in the formation of a 
strongly inclined inner starforming ring in the gSO galaxy 
having a radius of ~ 1 kpc. Its plane is orthogonal to the 
initial orbital plane of an interaction, i.e. its inclination to 
the gSO disc plane is about 60 deg. This ring is located in 
the inner bright part of the SO bulge and its contribution 
to the total stellar light (and, consequently, to the stellar 
kinematics derived from the spectra) in the central region of 
the galaxy is negligible (see left panel of Fig. I10|l . However, 
it is immediately revealed by the gas kinematics (see right 
panel of Fig. I10[) . as it is the only structure with significant 
gas content. Although the mass ratio of 1:10 in total mass 
is quite high, none of the retrograde encounters heated up 
significantly the large scale SO disc. 

The situation is totally different for prograde encoun- 
ters. They heat and thicken the gSO disc significantly, nearly 
destroying it for certain orbital types with high motion en- 
ergies. The gas and newly formed stars from the infalling 
satellite form a large scale co-rotating disc in the main plane 
of the gSO galaxy or between its initial plane orientation and 
the orbital plane if the large scale stellar disc is destroyed. 
If the large disc survives, the gas is concentrated at one of 
its resonances, probably the outer bar resonance, creating 
a starforming ring having R « 8 kpc quite thick (~ 1 kpc) 
in the radial direction (see Fig. Hip . As it resides in a re- 
gion where the surface density of the gSO disc is rather low, 
it creates significant contribution to the galaxy light. Since 
the initial metallicity of the interstellar medium in dSd is 
low, the newly formed stars in this ring will also have low 
metallicities being in agreement with young metal poor stars 
detected by our stellar population analysis in the outer ring 
of NGC 7217. Worth mentioning that due to the intense 
heating, the large scale gSO disc flares in the outer regions 
of a galaxy (see Fig |12|l , The outer disc must have formed by 
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Figure 11. Surface density of stars (grayscale) and gas (red) for 
a minor wet merger of gSO and dSd galaxies (mass ratio 10:1) on 
a prograde (direct) orbit. The orientation of the galaxy disc is 
chosen to match that of NGC 7217 (i = 30 deg). The axes are in 
kpc. 



this event not very long time ago, otherwise we would have 
expected to see young metal-poor stars formed in it migrat- 
ing inwards and erasing an observed quite sharp break in the 
age profile. The characteristic migration time of about 1 Gyr 
in a Milky Way-like galaxy using the mech anism proposed 
recent ly by iMinchev fc Famaevl (|2010l ) and iMinchev et al.l 
(|201lf ) allows us to constraint the formation epoch of the 
outer disc of NGC 7217. 

Hence, we conclude that the scenario of two consequent 
minor mergers at different orbits provides a plausible expla- 
nation to the observed structure of NGC 7217. 



6 SUMMARY 

We performed the analysis of internal kinematics and stellar 
population properties in the nearby isolated early-type disc 
galaxy NGC 7217 having star-forming rings by the NBuRSTS 
full spectral fitting of deep long-slit intermediate resolution 
spectroscopic data. 

The age and metallicity profiles of NGC 7217 exhibit 
three components that well correspo nd to those found in the 
struct ural analysis of the galaxy bv lSirchenko fe Afanasievl 
(2000): an old metal-rich bulge dominating at radii 0.4- 
1.2 kpc, the inner disc with the mean age of 5 Gyr and 
slightly subsolar metallicity at 1.6-4 kpc, and the outer 
metal-poor relatively young disc (t =2-3 Gyr) beyond 
4.8 kpc. 

We analysed the parametrized major- and minor-axis 
kinematical profiles of NGC 7217, recovered the velocity el- 
lipsoid in the regions dominated by the exponential discs 
and obtained the scaleheight radial profile. We concluded 
that the inner disc is thin (zq = 0.2 . . . 0.7 kpc) while the 
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Figure 10. Left panel: surface density of stars and gas (inset) for a minor wet merger of gSO and dSd galaxies (mass ratio 10:1) on a 
retrograde orbit. Right panel: line-of-sight velocities of stars and gas (inset) for the same merger. The signature of the inner polar ring 
seen edge-on is clearly visible. The orientation of the galaxy disc is chosen to match that of NGC 7217 (i = 30 deg). The axes are in kpc. 
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and showed that even the highest- density neutral-hydrogen 
ring at R = 6 kpc should be stable. 

From the analysis of emission line ratios in the resid- 
uals of the spectral fitting of NGC 7217, we ruled out a 
strong nuclear activity as the mechanism of circumnuclear 
gas excitation and concluded that the peculiar ionized-gas 
kinematics in the central region of NGC 7217 was due to the 
presence of the inner polar disc. 

By comparing our data with the numerical simulations 
of wet minor mergers in the GalMer database, we explained 
the observed structure and kinematics of NGC 7217 by the 
two consequent minor mergers with satellites having differ- 
ent initial orbits. The encounter on a retrograde orbit re- 
sulted in the formation of the inner polar disc, while the 
outer star-forming ring and the outer disc thickening were 
consistent with a minor merger on a prograde orbit. 
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Figure 12. Disc thickening by minor mergers. The edge-on view 
of the GalMer model of a giant SO galaxy in isolation is shown on 
top, the remnants of minor mergers with a gas-rich dSd galaxy 
(mass ratio 10:1, orbit type 1, see text) on retrograde and pro- 
grade orbits are displayed in the middle and in the bottom re- 
spectively. The contours are every 1 mag arcsec~ 2 , the axes are 
in kpc. 



outer disc flares reaching the half-thickness of 2.5 kpc in the 
region of the outer starforming ring. Using the two-fluid ap- 
proach, we evaluated the gaseous disc stability in NGC 7217 
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